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Abstract. The new generation of lithography tools use high energy EUV radiation 
which ionizes the present background gas due to photoionization. To predict and 
understand the long term impact on the highly delicate mirrors It is essential to 
characterize these kinds of EUV-induced plasmas. We measured the electron density 
evolution in argon gas during and just after irradiation by a short pulse of EUV light at 
13.5 nm by applying microwave cavity resonance spectroscopy. Dependencies on EUV 
pulse energy and gas pressure have been explored over a range relevant for industrial 
applications. 

Our experimental results show that the maximum reached electron density depends 
linearly on pulse energy. A quadratic dependence caused by photoionization and 
subsequent electron impact ionization by free electrons - is found from experiments 
where the gas pressure is varied. This is demonstrated by our theoretical estimates 
presented in this manuscript as well. 
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1. Introduction 

There is an increasing global demand for more computational power and memory 
capacity at lower cost and higher energy efficiency. The multibillion-dollar 
semiconductor industry is struggling to meet this demand and this is why they 
continuously aim to further miniaturize computer chips. To achieve this further 
miniaturization, the upcoming generation of commercial lithography tools use extreme 
ultraviolet (EUV) radiation at a wavelength of 13.5 nm, being equivalent to a photon 
energy of 92 eV. This EUV radiation will be generated by a pulsed light source based 
on plasmas in tin vapour, which are generated by irradiating liquid tin droplets with an 
intense laser beam BE]. The absorption cross section of EUV radiation in gasses is very 
high [3], so ideally the lithography tools would operate at high vacuum. However, for 
technical reasons manufacturers are forced to operate the tool with a low background 
gas (0.1-30Pa). Consequently, the interaction between the EUV photons and the 
background gas generates a plasma by photo-ionization, i.e. an EUV-induced plasma, 
everywhere it travels. In the vicinity of surfaces, the plasma induces strong electric 
helds due to the formation of a plasma sheath. Typically, the thickness of these sheaths 
is in the order of 5-10 mm and the electric helds are in the order of lO'^Vm”^ |1]. 
These strong electric helds accelerate ions to high velocities towards the highly delicate 
and expensive EUV multilayer optics, which can have a long term impact on their 
operation [5] . Eventually, controlling this long term impact requires more understanding 
of the fundamental processes in EUV-induced plasmas. 

Previous research used numerical simulations to investigate the fundamental 
properties of the plasma and its ehect on optical elements mu- In that study the 
authors also attempted to measure the electron density -being one of the key plasma 
parameters- using Langmuir probes; however, they concluded that the used probes are 
not feasible [6]. In ref. [7] we reported for the hrst time non-intrusive measurements 
of the electron density in an EUV-induced plasma in argon. For these measurements 
we applied a method known as microwave cavity resonance spectroscopy IMU. Where 
our previous work focused on solely the time evolution of the electron density and 
characterization of the MCRS method, the current paper explores the density over a 
range of parameters representable to those used in commercial lithography tools. 

The goal of this research is to answer the following question: How does the density 
of electrons in a plasma generated by pulsed EUV-radiation depend on gas pressure and 
the intensity of the EUV radiation? 

2. Microwave cavity resonance spectroscopy 

The applied diagnostics to measure the electron density is a method known as microwave 
cavity resonance spectroscopy (MCRS). This method is extensively used to measure 
the electron density in other types of plasmas HHII], while, recently, we published 
results of the hrst application of MCRS to EUV-induced plasmas [7]. Since the MCRS 
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technique is described in great detail in these publications, we will suffice with only 
briefly describing the key concepts here. In MCRS measurements a standing wave is 
excited in a (cylindrical) resonance cavity. This wave only exists at speciflc frequencies, 
called the resonant frequencies ujq, which depend, amongst others, on the permittivity 
of the medium inside the cavity. When free electrons are created in the cavity, the 
permittivity changes, due to which the resonant frequency will shift with Au. From 
this Acn the electron density can be determined [9]: 


He 


e2 uo 


( 1 ) 


with rUe the electron mass, e the elementary charge and ujq and oo the resonant 
frequencies without and with plasma, respectively. Note that this method gives the 
averaged electron density weighted with the square of the local electric held of the 
standing wave of the excited mode E{r) [9]: 
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To determine the local electron density, the spatial profile on the electron density needs 
to be known. Since this profile is unknown, we will present only the square-electric-fleld 
weighted averaged electron density. 


3. Experimental set-up 

We used a xenon-based EUV source to generate EUV radiation [1]. The experimental 
details of this source and the diagnostics will be discussed in this section. 

3.1. EUV source 

The pulsed EUV radiation, with a pulse duration of 150 ns and a repetition rate 
of 500 Hz, is generated by a pulsed xenon pinch plasma in the source chamber (see 
Figured]) [12]. The radiation is collected by a set of elliptic multilayer mirrors in the 
collector chamber, which focussed the EUV radiation in the intermediate focus (IF) in 
the measurement chamber. The EUV beam had a waist of 4 mm and a divergence of 
10°. To allow high vacuum in the measurement chamber, the collector and measurement 
chamber are differentially pumped and separated by an aperture of 4 mm. In front of the 
aperture, we placed a spectral purity Alter (SPF), which transmits between only 10 nm 
and 20 nm and, as such, prevents out-of-band radiation in the measurement chamber. 

3.2. Measurement vessel 

We placed diagnostics to study the induced plasma in the measurement vessel. As 
mentioned earlier, the main diagnostic technique was microwave cavity resonance 
spectroscopy. In addition to that, we determined the temporal shape of the EUV pulse 




Exploring the electron density in EUV-induced argon plasmas 


4 



Figure 1. (a) Schematic drawing of the experimental set-up. The EUV radiation is 
created in the source chamber by a xenon pinch discharge. The radiation is collected 
by a set of elliptic mirrors and focused in the measurement chamber. The cavity is 
placed in the focus. The gas environment of the measurement chamber is screened from 
the collector chamber by a spectral purify filter (SPF). This filter transmits between 
approximately 10-20nm. The argon pressure in the measurement chamber is varied 
with a needle valve, (b) Drawing of the microwave cavity with an indication of the 
EUV beam [7]. 


from the secondary electron emission from a copper disk. These techniques will be 
described in more detail below. 

We determined the spectrally integrated EUV power with a temperature sensor 
which heats up due to the absorption of EUV radiation. The temporal temperature 
increase is a measure for the time-averaged EUV power. The sensor consists of a copper 
disk (diameter of 25.4 mm and thickness of 1.0 mm) and a temperature transducer with 
a repeatability of ±0.1 °C. The temperature of the sensor is affected by the EUV power 
Pirn radiative losses and conductive losses. ^AT is then given by: 

IaT {AaAdT^ + X,)AT - QaA^T^ AT^ , (3) 

dt uicCp \ '■ -^^ ,—- / 

\ Cl Cl / 

where me is the mass of the copper disk, Cp is the specihe heat capacity of copper, a is 
the Stefan-Boltzmann constant, A^ is the area of the copper disk, Tq is the temperature 
of the environment and Ac is the conductive heat loss coefficient. During a power 
measurement, we hrst expose the power sensor to EUV radiation, due to which the 
temperature increases. After a few minutes of exposure, we close a shutter and the 
power sensor cools down. The resulting temperature curve is shown in Figure [2J We 
determined the heat loss terms by htting the temperature decrease with equation [3] 
{Pin = 0). With these -now known- loss terms we can £t the temperature increase and 
obtain the EUV power. The error in the power determination is 5 %. In the case of the 
sensor response of Figure [21 the power is (0.026 ± 0.001) W, which means an EUV pulse 
energy of (53 ± 3) pJ at a repetition rate of 500 Hz. 
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Figure 2. Heating and cooling of the temperature sensor (every 5th data point is 
shown to improve readability). The heat loss terms are determined from the cooling 
curve (Cl = (9.1 ± 0.1) x 10“^ WK“^ and C 2 = (3 ± 2) x 10“® WK“^). These terms 
are used to fit the heating curve with equation|3l The EUV power is (0.026 ± 0.001) W. 


To measure the temporal evolution of the radiation energy during an EUV pulse, 
we replaced the temperature sensor (copper disk and temperature transducer) by a 
copper disk only. This disk is connected to ground via a 50 U-resistor and we measured 
the voltage over this resistor with an oscilloscope. Due to the photo-electric effect, the 
EUV radiation induces a current. This is why, the temporal evolution of the current is 
representative for the temporal evolution of the EUV pulse. 

We varied the argon pressure in the measurement vessel between 0.1 Pa and 20 Pa 
with a needle valve. 

3.3. MCRS set-up 

We placed the microwave resonant cavity around the intermediate focus in the 
measurement vessel. This cavity is made of aluminium and has an inner radius of 33 mm 
and an inner height of 20 mm. It has entrance and exit holes for the EUV beam with 
radii of 6.5 mm. Inside the cavity two copper antennas are placed in the side-walls at 
opposite sides, one functioning as transmitter, the other as receiver. These antennas do 
not face the EUV light directly. A low power (10 mW) microwave generator (Stanford 
Research Systems SG386) excited the TMqio mode at (3.496 74 ± 0.000 04) GHz. We 
calculated the electric held of the TMqio mode in the exact cavity geometry with the 
plasimo platform [T^ITT] and the result is shown, along with the same resonant mode 
in a similar cavity without holes, in Figure |3l This graph shows that the entrance and 
exit hole have a small inhuence on the electric held prohle of the excited mode. The 
electric held is maximum in the centre of the cavity. This means that the electron 
density is most ehectively sampled in the centre of the cavity (see equation [2]), which 
is also the position where the plasma is created during and just after the EUV pulse 
when the narrow EUV beam is directed axially through the centre of the cavity. The 
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Figure 3. Modelled electric field of the TMgio mode in the experimental cavity [7]. 
The electric field is calculated with the plasimo platform [T^14) . Moveover, the electric 
field of an ideal cavity without holes is shown. 


response of the cavity is measured with a microwave detector (Hittiti HMC602LP4) 
which has a time response of 10 ns. The resonant frequency is determined by sweeping 
the microwave generator over a frequency range and measuring the response. At every 
point in time we £t the cavity response as function of the microwave frequency with 
a (two component) Fourier series. The frequency at which this £t is maximum, is the 
resonant frequency. It has been shown in that the response time of the cavity is 14 ns 
and that the combined response time of the cavity and the detector is about 17 ns. The 
error in the square-electric-field weighted averaged electron density is less than 30 % 
with a detection limit of 2 x 10^^ m“^, as has been shown as well in [7]. 

4. Results 
4-1- Gas pressure 

We measured the electron density as a function of time for various gas pressures in 
the range 0.5-15 Pa. Results are shown for respectively short and long time scales in 
Figures m and [5l These results show that the maximum weighted and averaged electron 
density increases with increasing pressure. At low pressures this maximum is reached 
at the end of the EUV pulse, of which the temporally resolved intensity is also plotted 
in Figure 0]), while for higher pressures it is reached later. Furthermore, the results 
indicate that the decay time becomes longer if the pressure increases. The noise level 
becomes signihcant below a density of 1 x 10^^ m“^ because we approach the detection 
limit of our experimental set-up ( rs-/ 10i2m 

The maximum electron density reached after the EUV pulse is plotted as a 
function of the gas pressure in Figure |6]for two different EUV pulse energies ((17 ± 1) pJ 
and (53±3)pJ). We could nicely £t the data with a quadratic function, as will be 
explained in section |5l For an EUV pulse energy of (53 ± 3) pJ the £t has the following 
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Figure 4. The electron density as a function of time for various gas pressures. The 
green dashed curve shows the normalized EUV intensity. The energy of the EUV pulse 
is (53 ± 3) pJ. 



time (ns) 

Figure 5. The electron density as a function of time for various gas pressures at a 
larger time scale. The energy of the EUV pulse is (53 ± 3) pJ. 


expression: 77 ,™“^ = (1.4 ± 0.5) x 10^^ + (1.5 ± 0.2) x 10^‘^p. The expression of the fit 
for an EUV pnlse energy of (17 ± 1) pJ is: = (2 ± 1) x 10^^p^ + (5.3 ± 0.9) x lO^^p. 

4-2. EUV pulse energy 

We also measnred the electron density as a fnnction of time for various EUV pulse 
energies. Results are shown in Figure [7] and demonstrate that the maximum electron 
density to be reached increases with increasing EUV pulse energy. Furthermore, they 
show that decay times do not signihcantly depend on the pulse energy. 

The maximum electron density reached is plotted as a function of the EUV pulse 
energy in Figure El We could fit the data well with a linear fit with an intercept fixed 
at zero. The slope of this £t is (2.1 ± 0.2) x 10^^ m“^ pj“^. 

























Exploring the electron density in EUV-induced argon plasmas 



Figure 6. Maximum electron density as a function of gas pressure for two different 
EUV energies, i.e. (53 ± 3) pJ and (17 ± 1) pJ. The expressions of the quadratic 
fits are, respectively, = (1.4 ±0.5) x lO^^p^ ± (1.5 ±0.2) x 10^'^p and 77 ,™““ = 

(2 ± 1) X lO^V ± (5.3 ± 0.9) X lO^V 



0 100 200 300 400 

time (|4S) 

Figure 7. The electron density as a function of time for various EUV pulse energies. 
The pressure was 5 Pa 


5. Discussion and interpretation 

We have measured the electron density as a function of gas pressure and EUV pulse 
energy. In this section we will discuss the results in therms of relevant physical processes. 

5.1. Gas pressure 

Our experimental results indicate a quadratic dependence of the maximum reached 
electron density on the gas pressure. To gain more insight we explore the production 
and loss processes of electrons in this conhguration. The most relevant among these are: 

(i) Production due to photo-ionization Ppi. 
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Figure 8. Maximum electron density as a function of EUV pulse energy. The fit has 
a slope of (2.1 ± 0.2) x and the intercept was fixed at zero. 

Table 1. Characteristic time scales for various processes in EUV induced plasmas. 


Process 

Typical time scale 

EUV pulse length 

100 ns 

Post EUV ionization time 

l~ 10 ps 

Time to reach wall 

10 -loops 

Total decay time 

0 . 01 - 1 ms 

Time between EUV pulses 

2 ms 


(ii) Production due to electron impact ionization Pgu 

(iii) Loss of initial high energy electrons L/e- 

(iv) Indirect loss due to averaging Lav 

(v) Loss due to recombination L^ec- 

Next, these processes will be described briefly. 

Photo-ionization The first production process is absorption of EUV photons 7 and 
subsequent photo-ionization Ppi j3]: 

7 -h Ar —> e~ -|- Ar^. (4) 

with a cross section at 92 eV of 1.1 x 10“^^ m^ |3] and an ionization energy of 
15.8 eV [ig or: 

7 -|- Ar —>• 2 e“ -I- Ar^^. (5) 

with a cross section at 92 eV of 2.5 x 10“^^ m^ [3] and an ionization energy of 
43.4 eV [15]. Triple ionization is negligible under our circumstances [3]. Due to 
momentum conservation, almost all excess energy is transferred to the electrons 
while the ions remain roughly at room temperature. So every absorbed photon 
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creates on average 1.2 electrons with an averaged energy of 66 eV (80% single 
ionization and 20% double ionization). 

Electron impact ionization Secondly we have production due to electron impact 
ionization Pei [iniin]: 


Ar —> 2 e + Ar"*" 
—)■ 3e“ + Ar^’*', 


( 6 ) 

(7) 


with cross sections at 66eV of 2.7 x 10 [16] and 1.3 x 10 [T7] . 

respectively. 

Loss of high energy electrons As mentioned before, the initial electrons have an 
average energy of 66 eV. These high-energy electrons escape from the centre of the 
cavity, leaving the positive ions behind. As a result a potential difference is induced 
between the plasma in the centre and the wall. The governing electric held conhnes 
the remaining electrons in the plasma. To estimate the charge density needed to 
conhne the plasma, we approximate the cavity with the initial plasma just after 
the EUV irradiation by the conhguration of a coaxial cable.The core of the coax 
represents the positive ion column while the shield of the coax represents the cavity 
wall. The electric held P in a coaxial cable can be determined with Gauss’s law: 

lE-dd=^, ( 8 ) 

Ja eo 

where Qenci is the charge enclosed by surface A. Solving this equation for E and 
integrating over r yields the potential between the core of the core and the shield: 


AG 


27rr£o \RcJ 


(9) 


with Xq the charge density per unit length, Rg the radius of the sheath and Rc the 
radius of the core. If we rewrite this equation for a plasma (A^ = Zen^rp, Rc = rp 
and Rs = R), the potential between the ions (core) and wall (shield) reads: 


AG = 


Zen+rp 

2^0 


In 



( 10 ) 


where Z is the ion charge, is the ion density in the core, is the radius of 
the plasma and R is the radius of the cavity. We assume that the radius of the 
core is equal to the radius of the EUV beam (2 mm). The average ion charge is 
1.2 (see previous point). From our calculation we hnd a necessary ion density of 
6 X 10^^ m“^ in order to conhne all electrons. To compare this to the measured 
electron density, we need to convert the ion density to a squared-electric-held 
weighted averaged ion density. Since we know the mode of the cavity (TMqio) 
and the local (relative) electric hels strength, we can convert the ion density to 
an averaged density using equation [2l which results in an averaged ion density of 
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6 X 10^^ m“^. Since the average ion charge is 1.2, this corresponds to an averaged 
electron density of 8 x 10 ^^ (initially quasi-neutral plasma), as to be measured 
with our MCRS method. So the loss term due to the fast electrons is small compared 
to the measured electron density and can be safely neglected. 

Indirect loss due to averaging Another loss process is related to the square-electric- 
held weighted averaging Lav If the plasma expands toward the wall, the apparent 
averaged density decreases since the E-£eld decreases, i.e. electrons are sampled to 
a lower extent. However, the E-£eld in our experimental cavity is approximately 
constant for 10 mm around the cavity axis (see Figure [3]). We estimate the time 
it takes the plasma to expand to 10 mm from the measurement at 0.5 Pa, since 
electron impact ionization is negligible at this pressure. The cavity-averaged density 
is constant for 1 ps (see Figure llj), so the loss due to averaging is not signihcant 
before 1 ps. 

Loss due to recombination Since it takes even longer for the plasma to reach the 
wall (radius of 33 mm) than the previously mentioned 1 ps, losses due to wall 
recombination L^ec are also neglected at time scales shorter than 1 ps. 


Typical time scales related to these processes are summarized in table [TJ 

The maximum electron density is always reached before 1 ps. Consequently, the 
rate equation for the electron density only includes Ppt and PeP 
dUe 
dt 


Ppi T Pei ^fe ^av 

'-V- 


=0 


1.2apinaEEuv(t) 


T hei‘^e'^a 


( 11 ) 


^EUVpphT 

= a{t) + bue, 

where Up, is the photo-ionization cross section, Ua is the background density, Aeuv is 
the cross section of the EUV beam, Eph is the photon energy, r is the duration of the 
EUV pulse, kei the electron impact ionization rate and 

Eeuv, 0 < t < r 
0, t > T 


^EUV 


(t) = 


( 12 ) 


where Eeuv is the EUV pulse energy. 
Appendix [71) : 


n, = 


^-2.(7piUaE EUV 


The solution of equation [TT] reads (see 


(13) 




AEuvEph 

This means that the maximum reached electron density increases linearly with EUV 
pulse energy and quadratically with the background density. This is exactly what we 
have observed in the experimental results in the framework of the pressure dependence 
study discussed in this section; we could £t a quadratic function perfectly through our 
experimental data (Figure [ 6 ]). 

The linear coefficient in the quadric hts in Figure [ 6 ] has a value of 
(1.5 ±0.2) X 10^"^m“^Pa“^ and (5.3 ± 0.9) x 10^^m“^Pa“^ for an EUV pulse energy 
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of 53 pJ and 17pJ, respectively. Note that, these values are valid for the squared- 
electric-held weighted averaged electron density. To convert these to coefficients which 
are valid for the peak electron density in the centre of the cavity, we need to assume a 
spatial distribution for the electron density just after the EUV pulse. We assumed that 
during the EUV pulse this distribution is a square with a width equal to the FWHM (full 
width at half maximum) of the EUV beam, i.e. 4 mm. Since we know the mode of the 
cavity (TMqio), we can convert the averaged electron density to a density in the centre 
of the cavity using equation |2l This results in linear coefficients of 1.4 x 10^® m“^ Pa“^ 
and 5 x 10^®m“^Pa“^ respectively. These coefficients correspond to the following term 
(linear coefficient) in equation [131 


at — 


l.^apiUaEEUV 1 
^EUvEph p’ 


which we can rewrite using the ideal gas law: 


at = 


1.2apiEEuv 

^EUvEphkbTg 


(14) 


(15) 


where p is the gas pressure, kb is the Boltzmann constant and Tg the gas temperature. 
We can estimate this linear coefficient from the absorbed EUV pulse energy. The cross 
section, photon energy and EUV energy depend on the wavelength of the photon A, so 
we write the previous equation as: 

^ ^■‘2o'pi{X)EEuvW 


at = 


-dX, 


( 16 ) 


kbTgAEuvhc/X 

where h is the Planck constant and c is the speed of light. To calculate this coefficient, we 
need the photo-ionization cross section and the EUV energy. The cross section is taken 
from Saito and Suzuki [3]. The spectrally resolved EUV pulse energy is calculated 
from the relative EUV spectrum recorded by Kieft [18] and a wavelength integrated 
energy measurement of the EUV source (53pJ per pulse). This results in a linear 
coefficient of 2 x 10^®m“^Pa“^, which corresponds well with the experimental value 
(1.4 X 10^® m“^Pa“^). When we decrease the EUV energy by a factor 3 to 17pJ, the 
theoretical coefficient also decreases by a factor 3 due to the linear dependence on the 
EUV pulse energy. This is also observed experimentally (a factor 2.8 ± 0.8). 

The difference in the quadratic part of the hts in Figure [6] when decreases the 
energy with a factor 3, should also be a factor of 3 since it scales with the EUV pulse 
energy (see equation [T^ . Experimentally we found a difference of a factor 7 ±4, which 
corresponds to 3 within the margin of error. 

Our experimental results (Figure Hj) show that at low pressures the maximum 
electron density is reached at the end of the EUV pulse, while at higher pressures 
this maximum is reached after well the EUV pulse. This difference in behaviour can be 
explained to the fact that at low pressures the total ionization is mainly governed by 
photo-ionization, hence the production of electrons stops at the end of the EUV pulse 
and the maximum is reached. At higher pressures electron impact ionization continues 
to create electrons after the EUV pulse at a signihcant rate. 
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5.1.1. Plasma decay As mentioned before, a potential difference is induced between the 
plasma and the cavity walls. On sufficiently long time scales, this causes an ambipolar 
flow ng. Due to this flow, the plasma expands to the walls of the cavity. Since the 
weighting factor of the MCRS diagnostic is smaller close to the walls of the cavity, 
the weighted average density decreases upon plasma expansion. However, it cannot be 
concluded that the total density also decreases. In general, the speed of this expansion 
is the ion acoustic speed {v = ^kBTe/mi) |20], which only holds after the hrst very fast 
electrons escaped to the wall. Since the electron temperature decreases with time and 
the rate of this decrease depends on the gas pressure, we can only indicate a range of 
expansion speeds: u = 5 x 10^ — 5 x 10^ms“^ for = O.l-lOeV. So it takes 7~70ps 
for the plasma to reach the wall of the cavity. When the plasma reaches the wall, the 
electrons and ions recombine. Typically, the decay time r of the plasma is then governed 
by the ambipolar flow IZH: 


T = 



(17) 


with R and H the radius and height of the cavity and Da the ambipolar diffusion 
coefficient. The equation above is only valid if an ambipolar equilibrium has been 
established. The ambipolar diffusion coefficient depends on the ion mobility /Xj, the 
electron temperature Tg and ion temperature Tj (in eV): 


Da — Pi ^Tg + I ( 18 ) 

The ion mobility decreases if the pressure increases |22l|23]. This perfectly explains our 
experimental observation of the plasma decaying faster at lower pressures. 


5.2. EUV pulse energy 

We also studied the dependence of the electron density on the EUV pulse energy at 
5 Pa. The experimental results (Figure E]) show that the maximum electron density 
depends linearly on the EUV pulse energy. The slope of squared-electric-held weighted 
averaged electron density against the pulse energy is (2.1 ± 0.2) x This 

linear dependence is also predicted by theory (equation [13]). 

The experiments (Figure [7|) show that the decay time does not depend on the EUV 
pulse energy. This is also clear from equations HT] and UHl 


6. Conclusions 

We experimentally explored the electron density evolution in an EUV induced plasma in 
argon using microwave cavity resonance spectroscopy. The dependence of the electron 
density on the EUV pulse energy and gas pressure was investigated. We found that 
the maximum reached electron density depends linearly on the EUV pulse energy as 
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also expected from theory. Also, our experimental results showed that the maximum 
reached electron density depends quadratically on the gas pressure. The linear term 
can be attributed to photo-ionization while the quadratic term is due to electron impact 
ionization. Finally, we found that the decay time of the EUV-induced plasma, where 
plasma recombination on the cavity walls is dominant, becomes longer if the pressure 
increases, which is explained by an increase of the ion mobility, hence the ambipolar 
diffusion coefficient, with increasing pressure. 

To gain more insight in the decay of the plasma, the spatial distribution of the 
plasma density as function of time should be determined. This can for instance be 
done by imaging the plasma. Furthermore, it is recommended to study the plasma in 
gasses with a smaller absorption cross section, such as hydrogen or helium, as used in 
lithography tools. 


7. Appendix 


The general solution of the equation: 

dUe _ l.2apinaEEuv{t) 
dt AEuvEphT 

= a{t) -f bue- 


+ keiUeU 


a 


is [24] : 

He = J 


where M{t) is the integrating factor: 


M{t) = 


/ 


bdt. 


(19) 

( 20 ) 

( 21 ) 


( 22 ) 


With the initial condition ne(0) = 0, this leads to the following solution: 


Heit) 


•)<*<- 


(23) 


where a* = a(0 < t < r). The maximum density in Figure 0] is reached at tm > r, so 
the equation for the maximum density is: 


^ ^ _ ^ 24 ) 

The term btm in the exponent is smaller than 1, so we can simplify the equation with a 
Taylor expansion: 

^max ^ ^ -F ^bHl^ - 1 - b{tra - t) - (25) 

(^1 + b{tm - ^r) 


= a T 


(26) 
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If we substitute a* and b we get the following expression for the maximum density: 


l.2a„inaE 


„ / / _ 

- 


pina-T^EUV 

AEuyEph 


1 -|- keiTta ( tfyi ^ T 


(27) 
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